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Abstract 



We show that the simultaneous determination of the leptonic CP-odd phase S and 
the angle ^13 from the subleading transitions Ve i^fi and De — > i^fi results generically, at 
^ • fixed neutrino energy and baseline, in two degenerate solutions. In light of this, we refine 

I I a previous analysis of the sensitivity to leptonic CP violation at a neutrino factory, in the 

LMA-MSW scenario, by exploring the full range of S and 613. Furthermore, we take into 
(-H I account the expected uncertainties on the solar and atmospheric oscillation parameters 

and in the average Earth matter density along the neutrino path. An intermediate base- 
line of 0(3000) km is still the best option to tackle CP violation, although a combination 
■ of two baselines turns out to be very important in resolving degeneracies. 



1 Introduction 

Solar parameters in the range of the large mixing ang le MSW solution (LMA-MSW) [|, |] 
are expected to affect sizeably the neutrino oscillation ||3| probabilities at terrestial distances. 
The discovery of CP-violation in the lepton sector might then be at reach 



Consider a neutrino factory from a muon storage ring ||T^, IT^ with muon energies of 



some dozens of GeV. In |11] a detailed study of the potential of such a neutrino factory in 
determining the leptonic CP-violating phase, 5, and the angle ^13 was performed, within the 
LMA-MSW. Just the range < 5 < 90° and ^13 > 1° was analyzed, though. It has been 



argued that the CP phase can be only determined up to a sign |17[. In order to clarify 
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this issue the full range —180° < 6 < 180° has to be explored. Besides, it is important to 
understand what is the lower value of ^13 at which the sensitivity to CP violation is lost. 

A second limitation of the analysis performed in jl^ is that the atmospheric parameters, 
023 and Am23, ^ ^^^^ ^ solar ones, 612 and Amfg, were assumed to be known. Neither 
an error estimate on the average Earth matter density was included. It has been recently 



pointed out |2C] that the errors in these parameters might modify drastically the conclusions 
reached in ITl 



For the atmospheric parameters, the main source of new information will come from 
proposed long baseline accelarator experiments (KEK, Minos, Opera) while for the solar 
ones lying in the LMA-MSW range, mostly from Kamland[p2[. Minos is expected to measure 



sin^ 2^23 and Am23 at the 10% level [^3|, in the range allowed by SuperKamiokande [ p^ . In 
the neutrino factory this uncertainty is expected to be improved to about 1% from muon 



disappearance measurements! 25, 26, 13]. Concerning the solar parameters in the LMA-MSW 
regime, sin^ 29i2 and Am^2> Kamland will be able to determine them to better than 10% for 



Amf2 > 10 ^ eV^ and sin^ 2^12 > 0.7 [27], well before the time of the neutrino factory. 

In this work we study, by means of an approximate analytical formula, the existence of 
degenerate solutions for the parameters (^13, 6): solutions that give the same probabilities for 



the fgf^ and I'e'^^i transitions than the one chosen by nature. We complete the analysis of ]11] 
by considering the full range of ^13 and 5, which reveals the existence of these degeneracies, 
and include the errors with which the atmospheric and solar parameters will be known at the 
time of the neutrino factory. Furthermore, we discuss the uncertainty on the average Earth 
electron density along the neutrino path, and study its relevance in the range 0(1-10%). 

In section 2, we analyse the sensitivity to CP violation in the full range —180° < 6 < 180° 
and clarify the correlation of 5 with the unknown parameter ^13. In section 3, the results of 
ref. ]11] on the simultaneous fits to ^13 and 5 are extended to the full range of 5 and ^13 < 1°. 
Section 4 describes the method to include in the analysis the errors on the atmospheric 
and solar parameters and on the Earth matter density, taking properly into account the 
corresponding correlations between neutrinos and antineutrinos, as well as between different 
energy bins and baselines (when more than one baseline is considered). The results of the 
fits including these errors and the discussion of the relative importance of each of them are 
then presented. Section 5 states our conclusions. 



2 Sensitivity to CP violation 

Given the high intensity expected at the neutrino factory, the effects of the solar mass differ- 
ence, Am^2) J^ot negligible over terrestrial distances in the LMA-MSW scenario, opening 
the way to observing leptonic CP violation. 

The best way to measure 5 and ^13 is through the subleading transitions i^e and 
They can be measured at a neutrino factory by searching for wrong-sign muons 



]l|, I while running in both polarities of the beam, i.e. and /x respectively. 

The exact oscillation probabilities in matter when no mass difference is neglected have 



been derived analytically in ]pq]. However, the physical implications of the formulae in ]28] are 
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Amf 



not easily derived. Defining Aij = ' ^ convenient and precise approximation is obtained 
by expanding to second order in the following small parameters: ^13, A12/A23, 1^12/ A. and 
A12 L. The result is (details of the calculation can be found in 
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where L is the baseline, Bzf = A13I and the matter parameter, A, is given in terms of 
the average electron number density, ne{L), as A = \^Gf ne{L), where the L-dependence 
will be taken from p9| . J is defined as 



J = cos ^13 sin 2^13 sin 2^23 sin 2^12. 



(2) 



In the limit ^ ^ 0, this expression reduces to the simple formulae in vacuum 
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In the following we will denote by atmospheric, P^}^, solar, P'^°\ and interference term, 



P^^p^, the three terms in eqs. (|T]). 
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Figure 1: Contours = P*"' (left) and PB*'^ = p^°^ (right) on the plane 6*13, Amf 2; /c"" 

the three reference baselines. Am23 = 3 x 10~^ eV^ and 612 = O23 = 45°. 

It is easy to show that 

implying two very different regimes. When ^13 is relatively large or Ami2 small, the probabil- 
ity is dominated by the atmospheric term, since -P"^^ ^ psoi ^ -yy^g ^-^j refer to this situation 
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as the atmospheric regime. Conversely, when ^13 is very smah or Am^g large, the solar term 
dominates P^°' ^ • This is the solar regime. Fig. |l| illustrates the separation between 
the two regimes on the plane Amf2 and ^13 for neutrinos and antineutrinos, as derived from 
eq. (0). The area to the right(left) of the curves corresponds to the atmospheric (solar) regime. 

We discuss next the subtleties in the measurent of CP violation in both regimes. 



2.1 Correlation between 6 and ^^13 

The oscillation probabilities of eq. (|l]), from whose measurement 6 could be extracted, depend 
as well on ^23) Am|3, ^12, Amf2, ^ and ^13. Uncertainties in the latter quantities can then 
hide the effect of CP violation. Although the first five of these parameters are expected to 
be known at the time of the neutrino factory with a good accuracy, ^13 might well remain 
unknown. It is essential then to understand whether the correlation between ^13 and 6 can 
be resolved in such a way that CP violation is measurable. The effect of the uncertainties in 
the remaining parameters will be analysed in section 5. 

Consider a single beam polarity and a fixed neutrino energy and baseline. The expansion 
of eq. (1^) to second order in ^13 leads to 

Pu.M^.p,) = X±dl, + y±^i3 cos i^±5 - ^) + , (5) 

with obvious assignations for the coefficients X and Y, which are independent of ^13 and 6. 
Note that the solar term P^"'' is the same for neutrinos and antineutrinos. 

Consider for instance Pi^^u^- The question is how many values of (^13, 5) give the same 
probability than some central values chosen by nature (^13 , 6) . 

This requirement can be solved simply for ^13 as a function of S: 

Y~ Ai3L 



^13 = -TT^COS 



2X^ V 2 
± Sqrt 



\2X+ V 2 JJ X 
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Eq. (P) is a curve| of equal probability on the plane (^13, 6), which for most of the parameter 
space spans the whole range of 6. It follows that, at any baseline, it is not possible to 
determine S from the measurement of wrong-sign muons at fixed neutrino energy with a 
single beam polarity. 

The analogous exercise for antineutrinos can be carried out resulting in a different equal- 
probability curve, with the substitutions in eq. (P): 6 — > —6, —>■ X_(yL). Assume 
now that both the neutrino and antineutrino oscillation probabilities have been measured, 
still at fixed (anti) neutrino energy and baseline. The question is if the two equal-probability 
curves intersect at values of {613,6) different from (^13,5). This condition implies equating 
eq. (P) to the corresponding one for antineutrinos and solving for S, for small ^13 > 0. The 
resulting equation is rather complicated, but simplifies considerably in the atmospheric and 
extreme solar regimes. 

^The sign has to be chosen so that 613 > 0. Discontinuities in S can arise only when the argument of the 
square root becomes negative. 
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2.1.1 Atmospheric regime 



In this regime it is safe to keep terms only up to first order in 1+ /X+ (YL / X_ ) in eq. @ . As 
a result only the solution of eq. (^) with + sign in front of the square root is acceptable since 
013 > 0. Eq. (6) simplifies to: 



h3 - ^"l3-^ 



Ai3L\ f- A13L 

cos — COS 



The equation for 6 is then obtained from equating eq. (0) for neutrinos to that for antineu- 
trinos. The problem ammounts to finding the roots of a function of 6 which is continuous 
and periodic. Since it must have at least one root at 5 = 5, by periodicity there must be at 
least a second root in the range —180° < 5 < 180°. 

The second solution for 5 in this approximation is: 

sin 5 — z cos 5 



sin 6 — sin 6 



cos — cos = 2z 7^ , (8) 



l + z2 

sin 6 — z cos 5 



where 

C+ A13L _ 1 y+ y_ , 
..^ta„^, ; C,.-(^±-). (9) 

The corresponding value of ^13 is: 

1 sin (5 - z cos ^ Ci - . A13L 



^13-^13-2 i + c: 

Only for the value of 6 satisfying 



sin^. (10) 



tan 6 = z (11) 

do the two solutions degenerate into one. Except for this particular point, there are two 
degenerate solutions with the penalty that, if nature has been perverse in her choice of 5, one 
solution may correspond to CP-conservation and its image not, and viceversa. 

In vacuum this is not the case. Eq. (^) in the vacuum limit: C_ — > or 2; ^ 00, gives 
6 = TT — 6 so that only for 6 = zb7r/2 there is no degeneracy. Then the two solutions either 
break or conserve CP. 

In Fig. (§), we show the result of 5 as a function of 6 for ^13 = 8° for three reference 
baselines together with the vacuum result. The difference between 6 and 6 is maximal close 
to 5 = 0°, 180°. 

It is interesting to consider the different impact of these degenerate solutions at different 
baselines. At short baselines, the oscillation probabilities for neutrinos and antineutrinos are 
approximately the same for two reasons: 1) the relative size of the sin 5 versus cos 5 term in 
eq. (|^) is tan(Ai3L/2) <^ 1, 2) matter effects are irrelevant with the solutions approaching 
the vacuum case |ll|, |l^. Indeed, the expansion of eq. (0) for A13L/2 <C 1 simplifies to 

ei3 ~ ^13 (cos 5- cos d) . (12) 
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Figure 2: Degenerate value of 5 as a function of true value 5, for Oi^ = 8° and three different 
baselines; the vacuum result 6 = it — 6 is also shown. 

The same equation holds for antineutrinos, since X^{Y+) = X-(Y-) in this approximation. 
The two equations have collapsed into one, and consequently we expect to find a continuum 
curve of solutions (^13, 5) of the approximate form given by eq. As the baseline increases 
the probabilities for neutrino and antineutrino oscillations start to differ, not only due to the 
term in sin 5, but also because of the matter effects. A shift in 6 cannot in general be then 
compensated in the neutrino and antineutrino probabilities by a common shift of ^13, and 
only the two-fold degeneracy discussed above survives. 



2.1.2 Solar regime 

In this regime the second term in eq. (|l|) dominates, although the first term cannot be 
neglected in the analysis of degenerate solutions even for very small values of ^13. The reason 
is that there exist, at fixed neutrino energy and baseline, a pair of values (613,6) at which 
the first and third terms in eq. (jl]) exactly compensate both for neutrinos and antineutrinos, 
in such a way that they are indistinguishable from the situation with ^13 = and any 6. It 
is easy to find these values by setting ^13 = in eq. (^) and in the equivalent equation for 
antineutrinos. 6 is the solution of: 

tan 6 = — , (13) 

z 

and the corresponding ^13 is: 
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Taking as an example Amj^ = 3 x lO^^ eV^, L = 2810 km and = 0.3E^,E^ = 50 GeV, 
this point is: 

%~1.5°, (5~-165°. (15) 

As we will see in the next section, this solution is clearly seen in the fits. Alike to the pattern 
in the atmospheric regime, this degeneracy occurs only at fixed neutrino energy and baseline. 

In summary, this analysis implies that, even with the information from both beam polari- 
ties, there are in general two equally probable solutions, at fixed neutrino energy and baseline, 
for the parameters ^13 and 5. This is of course applicable to any experiment that tries to 
measure simultaneously ^13 and 6 from the subleading transitions Uet^^i and ^'e^'^i- These de- 
generacies can in principle be resolved by exploiting the energy and baseline dependence of 
the oscillation signals, as will be discussed later on. We point out that a supplementary mea- 
surement of the channel t'ei^rl'^ei^r) could also be of great help in resolving the degeneracies, 
if performed with enough accuracy, as the (5-dependent terms in the oscillation probabilities 
have the opposite sign than for the i^e^nii^e^ii) transitions. The experimental challenge of 
measuring those transitions is however much greater and we will not follow this avenue here. 



3 Simultaneous determination of 5 and 9 
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In 1 11] a detailed study of the possibility to determine simultaneously ^13 and 5 was performed, 
but only the range ^13 > 1° was explored for the LMA-MSW solution. According to Fig. ||, 
this corresponds to the atmospheric regime. However, the sensitivity to ^13 for the SMA- 
MSW solution (where only the term P"^*"^ survives) was found to extend to smaller values: 
013 ~ 0.15° at Am|3 = 3 x 10"^ eV^. For the LMA-MSW scenario, such a low value of 6I13 
lies well inside what we have named the solar regime, which is thus also relevant in the study 
of CP-violation for 6*13 < 1°. 

It is very illustrative to consider the different behaviour of the CP asymmetries, defined 
in [^0| , |5|, 0, in both regimes, as an indication of the sensitivity to CP violation. 

The (5-dependent terms in the oscillation probabilities of eq. (Q) are linearly suppressed 
on the two small parameters: Amf2 and ^13. In the atmospheric regime the leading term, 
patm ^ does not depend on Amf2, while it is quadratically dependent on ^13. The sensitivity 
to CP violation decreases thus linearly with Am^g while it is rather stable as ^13 decreases ^. 
There is of course a limit to this CP-insensitivity to ^13 and this is precisely when we enter 
the solar regime. In the solar regime the role of both parameters is interchanged and while 
the sensitivity to CP violation decreases linearly with ^13, it remains rather flat with Am^2- 

This behaviour is shown in Fig. |^, which displays the significance of the CP-odd asymme- 
tries defined in [ |30| , ^, ^ as a function of ^13, for Amf2 fixed, and viceversa. In both figures 
the change in behaviour coincides roughly with the limit between the atmospheric and solar 
regimes. These asymmetries have been obtained for a muon beam of 50 GeV providing 10^^ 
useful /i"*" and ii~ decays and a 40 Kton magnetized iron detector [^ ], which is our working 
setup in the present work, as it was in |ll 



^Although the relative importance of the terms in 5 with respect to the leading term increases with de- 
creasing fel, the statistical significance remains constant [M. 
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Figure 3: Significance of the asymmetries at L = 2810 km as a function of 613 for E"^ = 50 
GeV, Aml2 = eV^ (left figure) and Amj2 for 613 = 0.5° (right figure) . Am^g = 3x10"^ 
eV^ and 612 = G23 = 45°. The five dashed curves correspond to neutrino energy bins of 10 
GeV, and the solid curve is the average. 

In |jll|] (see also [jl^) 5 and ^13 were extracted from a spectral fit of wrong-sign muons 
signals for both polarities, using the exact oscillation probabilities in the approximation of 
constant Earth matter density. Note that an analysis based on CP asymmetries alone does 
not use all the available information. The observables used were the number of wrong-sign 
muons in five bins of energy for both beam polarities: 

A^±, (16) 

where A labels the baseline, i the energy bin and it the sign of the decaying muons. These 
numbers are given by: 

i-Ei+AE 

K^,±= $,(p)(E„L) a,(p)(E,) P,(p)(£;„L,0i3,<5,a) (17) 

J Ei 

where a is the set of remaining oscillation parameters: ^23; ^12) Am^g, Amf2 and the matter 
parameter A, which were taken as known. denote the neutrino fluxes and the DIS 

cross sections. 

Simultaneous fits of the parameters 6 and ^13 were performed for three reference 
baselines L = 732 km (the Cern to Gran Sasso distance or Fermilab to Soudan), 3500 km (to 
be replaced in this work by 2810 km, the distance from Cern to La Palma) and 7332 km, as 
well as for various combinations of them. 

The cit 9- fixed baseline is of the generic form: 

= E E «P - ^Kp)Cii - ^^,P') ' (18) 

P,P' 

where C is the 2Ni,in x 2Ni,in covariance matrix, n^^ are the simulated "data" obtained from 
a Gaussian or Poisson smearing including backgrounds and efficiencies (for more details we 
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refer to [^)- In the combination of two baselines we have: 

xlx' = E E E «P - Kp)Ci:-,p-.i',,,pHp' - Kv) ' (19) 

1,1' i,j p,p' 

where C is now a matrix of dimension 4:N{jin ^ ^^Un 

As the errors on the a were neglected, C contained only statistical errors, 5n\p, which 
were assumed to be independent for different i, p and /: 

Cl,i,p;l',j,p' = ^ij ^pp' (f^^lp)^- (20) 

In |11| only the restricted range < (5 < 90° and ^13 > 1° was explored. For this 
reason the degeneracies appearing when the full range of 5 is considered were missed. We 
have repeated the analysis by considering the full range of 5 and ^13. In all the plots that 
follow realistic efficiencies and backgrounds have been included. We have checked that the 
degenerate images in the plots shown below appear at the points indicated by the analysis 
of the previous section. The false images are somewhat softened, as different neutrino and 
antineutrino energies enter in the analysis, but still visible. We present fits only for Am23 > 0. 
The opposite case gives better results: for Am23 < 0, the statistics for the signals of positive 
and negative wrong-sign muons are closer (so that the difference is more neatly seen) because 
matter effects enhance in this case i^p^p^, compensating to a large extent the difference in the 
neutrino DIS cross sections Uy ~ 2(Tp. 

All the results shown below correspond to central values of the parameters in the LMA- 
MSW scenario: Amfg = lO""^ eV^, /^ml^ = 3 x 10"^ eV^ and 6*12 = 6*23 = 45°, except in Fig. 



12, where the full range of Amf2 is considered 



3.1 Atmospheric regime 

Recall that the pure CP-effects peak at distances in the range 2000-4000 km [^, |^ |ll[. In 
Figs. Q we show the results of the fits including efficiencies and backgrounds for L = 2810 
km for central values of ^ = -90°,0°,90°, 180° and for ^13 = 2° (left) and ^13 = 8° (right). 
The energy dependence of the signals is not significant enough (with our setup) to resolve 
the expected two-fold degeneracy. The second solution is clearly seen for the central value of 
(5 = 0° as an isolated island. For the central values 6 = —90° and 6 = 90°, the degeneracy 
is responsible for the rather large contours which encompass the two solutions. Notice that 
as ^13 diminishes the fake solution for 5 = 90° moves towards 6 = 180°, as it should (recall 
that, in the solar regime, the vacuum fake image lies at 6 = 180°). 

Figs. ^ show the fits for ^13 = 8° at L = 732 km and 7332 km. In the former, the expected 



continous line of solutions of the form given by eq. (12) is clearly seen. The measurement of 
6 is thus impossible at this baseline if ^13 is unknown. In the longer baseline, the sensitivity 
to 5 is similarly lost but for a different reason: the CP-signal is fading away (indeed the 
underlying degenerate solutions become much closer in ^13) and statistics is diminishing. 

In Fig. ^ we show the result of combining any two baselines. The two- fold degeneracy does 
not disappear completely in the combination of the two shorter baselines, while it does in the 
remaining two combinations. It is interesting that, while the shorter and longer baselines by 
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Figure 4: Simultaneous fits of 5 and 613 at L = 2810 km for different central values (indicated 
by the stars) of 6 = -90°,0°,90°, 180° and 613 = 2° (left), 8° (right). The value of 6 for the 
degenerate solutions is also indicated. 

themselves are not appropiate for CP studies, their combination is quite promising, due to 
the very different pattern shown by the correlation between ^13 and 5 in each of them. The 
overall conclusion is that the combination of any two distances is interesting regarding CP 
violation. 



3.2 Solar regime 

In Fig. we show the results of the fits including efficiencies and backgrounds for L = 2810 km 
for central values of 6 = -90°, 0°, 90°, 180° and for ^13 = 0.3° (left) and ^13 = 0.6° (right). 
Consider for instance the case ^13 = 0.3° (left): the degenerate images of the four points 
chosen appear grouped at the right/lower side of the figure. These are the solutions that 



mimic ^13 = as predicted in section |2.1.2 . The position of these solutions can be accurately 
predicted from the analysis of the previous section. The comparison of these figures with 
Fig. 1^ illustrates the expected decrease of the sensitivity to CP violation for very small ^13, 
consistent with the behaviour of the significance of the CP-asymmetry in the solar regime 
shown in Fig. ^. Note that at ^13 = 0.3, the sensitivity to CP is already lost for 5 = —90°. 

In Fig. P we show the result of combining the baselines L = 2810 km and L = 732 km 
with the longer one. As in the atmospheric regime, the degeneracies are nicely resolved in the 
combination of the intermediate and long baselines. However, in this case the combination 
of the short and long baselines is not as good. Note the degenerate solution that survives 
for 6 = 90° which appears centered at 5 = —90°. This is the result of the intersection of the 
approximately vertical contours found at L = 7332 km with the cos 5 contours at L = 732 
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Figure 5: Simultaneous fits of 6 and 613 at L = 732 km (left) and L = 7332 km (right) for 
different central values of 6 = -90°, 0°, 90°, 180° and 613 = 8°. 

km (see Fig. 0). 



New analysis: inclusion of expected errors on oscillation 
parameters and matter density 



Consider the uncertainty in the theoretical parameters a. These errors induce in general a 
correlation between the obse 
C of eq. (p^ is of the form: 



correlation between the observables, N^"^, which has to be taken into account. The matrix 



C^,p-v,j,p' ^ ^iv V ('^^Ip)' + E ^^^^ 

a 

where a{a) is the la error on the parameter a. 

Recent analysis of the expected uncertainty in the knowledge of the atmospheric param- 
eters at the neutrino factory indicate a ~ 1% uncertainty in Am23 and sin^ 2^23 1 25, |2^0 For 



the solar parameters we include the results of the analyses of the Kamland reach [^7|: 2% 
error in Amf2 and ±0.04 in sin^20i2, for maximal ^12, both at \o. Slightly smaller errors 
could be obtained for the product (sin 2^12 Amf2)^ which is the combination entering the 
relevant oscillation probabilities |^^, a refinement we neglect here. For the uncertainty on 
the matter parameter, we could not find any estimate in the literature. The dispersion 



Although these analyses have been done for the SMA-MSW solution or assuming that the solar parameters 
are known, we will assume that in the LMA-MSW scenario the errors on the solar parameters or in the matter 
term do not change this result. 
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Figure 6: Fits of b and ^13 combining any two baselines. 



of the different models of the Earth density profile |33] indicates an uncertainty of 1-2% for 
trajectories which do not cross the core, though. We consider a range between 1-10% for 
illustration. 
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Figure 7: Simultaneous fits of 5 and 613 at L = 2810 km, for different central values of 
S = -90°, 0°, 90°, 180° and ^13 = 0.3° (left), 0.6° (right). The value of 6 for the degenerate 
solutions is indicated. 
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Figure 8: Simultaneous fits of 6 and 613 at two combinations of baselines, L = 2810 + 7332 
km (left) and L = 732 + 7332 km (right), for ^13 = 0.6° and different central values of 
5 = -90°, 0°, 90°, 180°. 



The most important changes result from the uncertainty in 623 and in the matter param- 
eter A, with the former affecting mainly the measurement of ^13 and the latter the sensitivity 



13 



to 5. 



Recall Fig. 4, where no errors on the oscillation and matter parameters were included. 
Fig. ^ (left) depicts the results for 6 = 90° and —90° at L = 2810 km, including all errors (with 
an error in the matter parameter of 1%) compared (right) with the situation in which only 
the error on the atmospheric angle 623 is included. The two graphics in this figure are almost 
identical, showing that the dominant error is that of ^23- It affects mainly the determination 
of 013, a fact easy to understand: while the measurement of the leading transition i^^ 
is sensitive to sin^2023i and this can be measured with a 1% uncertainty, the subleading 
transition Vf, is proportional to sin^ 622,- For maximal ^23 mixing a 1% relative error 

in the former translates into a 6% in the latter. This is then the largest relative error of all 
the parameters that enter P'^*™^ which is dominant in most of the parameter space. Note 
that the effect of the errors is more important for larger ^13. In Fig. IC we show the results 
for the best combination of baselines when all errors have been included. The resolution of 
the degeneracies discussed in the previous sections is still achieved, but the contours have 
become sizeably larger. 
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Figure 9: Fits of 6 and 613 for various central values of 5 and 613 at L = 2810 km including 
all the errors on the remaining parameters (left plot) with AA/A = 1% and including only 
the error on 623 (right plot). 

The uncertainty in A is more relevant for the determination of 6, although if this error is 
controlled at the percent level the effect is negligible. In general, even for AA/A = 10%, the 
effect is far less important than the error induced by correlations between ^13 and 5. This can 
be seen with just one baseline, where the degeneracies survive. As an illustration, in Fig. ^ 
for L=2810 km, AA/A is varied in the range 1 — 10%, with 6 = 90°, ^13 = 8° and the errors 
on the remaining oscillation parameters included. This is to be compared with Fig. Q (right) 
where no uncertainties were assumed. The error in 6 is seen to be mostly dominated by the 
correlation of ^13 and 6. For the combination L = 2810 + 7332 km, where the degeneracies 
are resolved, the effect of the error in A is more important in relative terms (compare Fig. ^ 
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Figure 10: Fits of 6 and ^13 combining the two baselines: 2810 km and 7332 km, for various 
central values of 6 and ^13 including all the errors on the remaining parameters with /S.A/A = 



70. 



with the middle plot in Fig. P): a 10% error in A is a 50% increase of the error in 5. It is 
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Figure 11: Results from a simultaneous fits to ^13, 5 for 6 = 90° including errors on the 
oscillation parameters and varying the error on the average Earth matter density from 1% = 
10%. 

clearly desirable to have a per cent control over the average Earth matter density, which does 
not look a priori unrealistic. 

Recently, the authors of have presented an analysis of the sensitivity to CP violation, 
including the errors on the oscillation and matter parameters, with quite different conclusions. 
In particular, they state that CP violation can only be measured in a small window at the 
shorter baselines. There are a number of differences between their analysis and ours. They do 
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not consider the simultaneous determination of ^13 and 5, but include and ad hoc 10% error 
on ^13. They assume as well larger errors on the remaining solar and atmospheric parameters: 
a democratic 10%. We have seen that estimates from Kamland and the expectation from 
disappearance measurements at the neutrino factory give more optimistic results. Finally, 
they do not take into account correlations of the errors that these parameters induce on 
the different observables: wrong-sign muons in different energy bins and different polarities. 
Neither they include experimental background and efficiencies. 

Let us turn now to the case of smaller values of Amf2 allowed in the LMA-MSW range. 
For fixed ^13 , the sensitivity to 5 decreases linearly with Amf 2 ™ atmospheric regime and 
more slowly in the solar one. For the plots in the atmospheric regime of this section it does 
not necessarily imply, though, a linear scaling (with Am^g) of the error in 5, as the latter 
is mostly dominated by the existence of degenerate solutions on the plane (^13, (5), whose 
separation dos not follow such a linear pattern. 

It is interesting to understand how much of the LMA-MSW range can be covered in the 



discovery of CP violation. This is illustrated in Fig. |1J with a rough exclusion plot. For the 
hypothetical nature values 6 = 90° and the best combination of baselines, L = 2810+7332 km, 
the line corresponds to the minimum value of Amf2 at which the 99%CL error on the phase 
reaches 90° degrees, and is thus indistinguishable from 0° or 180° (i.e. no CP violation). The 
error on the phase is computed by taking the longest vertical size (upwards or downwards, 
whichever is longest) of the 99%CL contour from 90°. All errors on the parameters have been 
included. With this definition, there is sensitivity to CP violation for ^13 > few tenths of 
degree and Amf2 > 3 x 10"^ eV^ 

An analogous plot in [|l^] indicated better sensitivity. Several reasons account for the 
difference. First, only the possibility of separating 90° from 0° (instead of 180°, which is 
more constraining |20|) was considered there. Second, the correlation between ^13 and 5 was 



not taken into account: in order words, ^13 was fixed at its true value ^13 and only then the 
error on 6 considered. Finally, the errors on the remaining oscillation parameters and A were 
not included, although this makes a very small difference if AA/A = 1%. 

A final comment. As stated before, the problem of the correlations between ^13 and 6 
has to be faced by any experiment measuring just the and i>e ^ i^^i transitions. In 

particular this applies to the so-called "superbeams" |34]: intense neutrino beams from pion 
(and kaon) decay. They could provide, though, very useful complementary information to 



the neutrino factory in disentangling ^13 and 6, for their expected reach ^13 > 3° [34]. 



5 Conclusions 

A neutrino factory from muon storage rings, with muon energies of a few dozen GeV, is 
an appropiate facility to discover leptonic CP violation through wrong-sign muon searches. 
This requires that the solution to the neutrino solar deficit is confirmed to lie in the LMA- 
MSW regime, and the angle ^13 is larger than a few tenths of degree. Within this range, 
the sensitivity to CP-violation is lost only for the smaller values of the solar mass difference 
allowed by the LMA-MSW scenario. 

At the hypothetical time of the neutrino factory, the value of the parameters ^13 and 6 
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Figure 12: Sensitivity reach for CP violation as defined in the text on the plane (Amf2j^i3) 
for the combination of baselines L = 2810 and 7332 km. All errors are included. 

may be still unknown and will have to be simultaneously measured. In this paper we have 
considered the full range of possible values of 5. A relevant problem unearthed is the generic 
existence, at a given (anti)neutrino energy and fixed baseline, of a second value of the set 
{013,6) which gives the same oscillation probabilities for neutrinos and antineutrinos than 
the true value chosen by nature. It is a generic challenge for any future facility. The spectral 
analysis and the combination of baselines satisfactorily resolves this degeneracy. 

Furthermore, we have included in the analysis the expected uncertainty on the knowledge 
of the rest of the oscillation parameters (sin^023, Am^a, sin^Ou, Amf2) and on the Earth 
electron density. Noticeable changes result from the error on ^235 which affects mainly the 
uncertainty in ^13, and from the uncertainty on the Earth matter profile, which affects mainly 
the extraction of 5. The latter uncertainty is of little consequence if at the level of a few 
percent. It seems pertinent to us that a detailed geological analysis of the planned baselines 
is performed, to reinforce the expectations as regards CP violation. 

Realistic background detection errors and experimental efficiencies have been included in 
the analysis. The overall conclusion is that the optimal distance for studying CP-violation 
effects with neutrino energies of few dozens of GeV is still of 0(3000) km, although the 
combination of two baselines, one of which being preferably a very long one, is very important 
in resolving degeneracies. 
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